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METHOD FOR CREATING POLYNUCLEOTIDE AND 
POLYPEPTIDE SEQUENCES 

5 

CROSS-REFERENCES TO RELATED APPLICATIONS 
This application derives priority from USSN 60/067908, filed 
December 8, 1997, which is incorporated by reference in its entirety for all purposes. 

10 TECHNICAL FIELD 

The invention resides in the technical field of genetics, and more 
specifically, forced molecular evolution of polynucleotides to acquire desired properties. 



BACKGROUND 

15 A variety of approaches, including rational design and directed evolution, 

have been used to optimize protein functions (1,2). The choice of approach for a given 
optimization problem depends, in part, on the degree of understanding of the relationships 
between sequence, structure and function. Rational redesign typically requires extensive 
knowledge of a structure- function relationship. Directed evolution requires little or no 
20 specific knowledge about structure-function relationship; rather, the essential features is a 
means to evaluate the function to be optimized. Directed evolution involves the 
generation of libraries of mutant molecules followed by selection or screening for the 
desired function. Gene products which show improvement with respect to the desired 
property or set of properties are identified by selection or screening. The gene(s) 
25 encoding those products can be subjected to further cycles of the process in order to 

accumulate beneficial mutations. This evolution can involve few or many generations, 
depending on how far one wishes to progress and the effects of mutations typically 
observed in each generation. Such approaches have been used to create novel functional 
nucleic acids (3, 4), peptides and other small molecules (3), antibodies (3), as well as 
30 enzymes and other proteins (5, 6, 7). These procedures are fairly tolerant to inaccuracies 
and noise in the function evaluation (7). 
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Several publications have discussed the role of gene recombination in 
directed evolution (see WO 97/07205. WO 98/42727, US 5807723, US 5,721367, 
US 5,776.744 and WO 98/41645 US 5.81 1,238. WO 98/41622, WO 98/41623. and 
US 5.093.257). 

5 A PCR-based group of recombination methods consists of DNA shuffling 

[5, 6], staggered extension process [89, 90] and random-priming recombination [87]. 
Such methods typically involve synthesis of significant amounts of DNA during 
assembly/recombination step and subsequent amplification of the final products and the 
efficiency of amplification decreases with gene size increase. 
10 Yeast cells, which possess an active system for homologous 

recombination, have been used for in vivo recombination. Cells transformed with a 
vector and partially overlapping inserts efficiently join the inserts together in the regions 
of homology and restore a functional, covalently-closed plasmid [91]. This method does 
not require PCR amplification at any stage of recombination and therefore is free from the 
15 size considerations inherent in this method. However, the number of crossovers 

introduced in one recombination event is limited by the efficiency of transformation of 
one cell with multiple inserts. Other in vivo recombination methods entail recombination 
between two parental genes cloned on the same plasmid in a tandem orientation. One 
method relies on homologous recombination machinery of bacterial cells to produce 
20 chimeric genes [92], A first gene in the tandem provides the N-terminal part of the target 
protein, and a second provides the C-terminal part. However, only one crossover can be 
generated by this approach. Another in vivo recombination method uses the same tandem 
organization of substrates in a vector [93]. Before transformation into E. coli cells, 
plasmids are linearized by endonuclease digestion between the parental sequences. 
25 Recombination is performed in vivo by the enzymes responsible for double-strand break 
repair. The ends of linear molecules are degraded by a 5*->T exonuclease activity, 
followed by annealing of complementary single-strand 3' ends and restoration of the 
double-strand plasmid [94]. This method has similar advantages and disadvantages of 
tandem recombination on circular plasmid. 

30 

SUMMARY OF THE INVENTION 
The invention provides methods for evolving a polynucleotide toward 
acquisition of a desired property. Such methods entail incubating a population of 
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a .^anions to eenerate annealed polynucleotides 
compnsesheteroduplexes. lnen noWnuc leotide vanants or 

recombinedpolynucleot.de vanants. The resulting p . 

hp nrenared in the form of cellular extracts. 

introduced into host cells. The heteroduplexes are thus exposed to the host cells DN A 
repai r system in vivo. ^ ^ ^ ^ ^ 

c an he achieved, for example, by prov.dmga „ . comp onen, of a 

of a firs, vector, and a second poiynuc.eot.de vanan, s P„ d ^ 
5 sec ondvec,o, Thefirstandsecondvectorsarecon 

th e first and second polynucleotide variants occur ^ ^ ,„ ^ 

lineaI firs, vector, single-stranded fonns ^ foms of the fo , 

each other to form linear second vector, *td s,n ^ ^ . n 

20 and second vectors a™*, with each to form a - ^ * het£rodup , exes 

25 ^ reS ^t "thods, polynucleotide variants are provided in double stranded 
in some m j rf f„ m before the annealing step. Opttonally, 

form and are converted to smgle stranded form before 
such conversion.sby conducting 

^ v,u branded form, and the method further 
second polynucleotides is provided ,n double stranded 



PCT/US9 8/25698 



WO 99/39902 



10 



15 



20 



• rfcP first and second polynucleotides as components of first and 
linearized vector reannealwitneacnomci 

hneanzea vcu hnear s£Cond 

forms of the second linked vector reanneal - * ^ 
vector, and single-stranded linearized « ^ ^ 

step selects for transformed cells comprises m 

u a th P nolvnucleotide variants encode vanants of a 
„ m ,. nNA in some methods, the poiynucicouu 
chromosomal DN A. In so ^ . nolvnuc leotide variants comprises at 

p^-7j^«S^^— ----- 

least 20 vanants. In some metnoas, mc y f 

PCR or cassette mutagenesis. In some metnou*, 

, „ bacterial cells In some methods, the population of variant 
are mtroduced are bacterial cell, Qtldes having at lea st 90% sequence 

polynucleotide variants comprises at least i poiyn 

ideM " y ^ Tol — runner comprise a step of atleas, partially demethylating 
variant polynucleotides. Demelhylation can be achieved by PCR amplification 
passaging variants through methylauon-deficient host cells. 

' Some methods include a further step of sealing one or more nicKs ,„ 

h eteroduple* molecules before .posing the heteroses DMA repatr system. 

u j Trinucleotide variant is screened to 

lherrty * ,n some methods, the recombinant protein or secondary metabolite ,s 

formulated with a earner to form a pharmaceutical composuion. 
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oxidases, peroxidases and phytases. In 

a poiypeptide seiected from lh= group ^"^^ ormones, 
5 ^statin, —in, thymosin, -^rirn cJldoZn, ***** 

, 0 p o ly „ne,eo t ide S eneodesap^c f en, m es 

BRIEF DESCRIPTION OF THE DRAWINGS 
FiBure . iUustrates the process of he.erodup.ex formation using 

„ <PCR) wi,h one sc. of primers for each different sequence to 
polymerase chain reaction (PCR) witn on 

enzymes to Hne*nz= the target -I— formation usin g asymmetric 

Fig ure 3 iUustra.es a process of he«rod ^ ^ 

« singie primer poiymerase chain reaction (PCR) with 

-^'l^rrZ — of^que— nshes 
u , ilte ed to lineage pMJO and pRM40. ^ ^ 

flaA genes containing sequence information from both parent genes. 
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Fig ^ 9 sho»s ft. physica. maps of Aa^es u«*~* KB 
„ — wiin enhanced specrf.c aciiviry «■> » PM7-2 for Muran, 7* - 0» 

5 mu „ions ,„ — 7-2 and « » = CB » ^ 

Figure llshows posiuonsofDNAbasechangesandannnoacd 
,. s,— ns in — ECB deacyiase Muran, .5 — » -*«~ 

^a^ne^M-Usenesc.a.ed^.v,,.— e—onfoiioW 
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by in vivo repair (SEQ ID NO:7). 

Figure 14 illustrates the process used for Example 

• p r 1 and RC2 to yield thermostable subtilisin E. 
mutations m RC1 and RCZ y dones 

Figure 15 illustrates the sequences of RC1 and K^z 

v. fnm^nts of the reaction products of duplex formation as 
plck ed random* ^ ^^t, base posLns - ^ «=' - 

^07 clspohds « an amino acid subsutuuon ai 2, 8 in Ok ***** P™,-~ 
Figure ,6sh O w S *= r =su taO fs™400c 1 on esft o m .he W 

„ea,ed by heterodluplex formation and repair for rn.ua, ac.iv.ry <*> "» 
I Tn r a,,o Ai/ A, «. used ,o esrimare the enzymes' -""^^^ 
c L varianis are sorted and piotted in descending order. App—y 2.9 ,f the 
active vm double mutant containing both the 

clones exhibit a phenotype corresponding to the double m 

N181D and the N218S mutations. 

DEFINITIONS 

Screening ,s. in generai, a two-step process .n which one first physical 

property Selection is „ lection marker, which, in some genetic 

are achieved simultaneously bv expression of a selection mar 
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^ marker to survive while other cells die (or vice 
circumstances, allows cells expressing the marker to surv 

inrinrfe luciferase. £galactosidase and green 
versa> Exemplarvscreeningmemoersincludeluciierase.ps 

n: L , pin. Se.ection marxers include drug and -in resistance genes. ^ 
^taneous selection can an, do. ..cur i» the course of natural evolution. ,n th = present 
methods selection is performed by man. 

An exogenous DNA segment ,s one foreign (or heterologous, to the ceU or 
homologous to the cell hut in a position within the host ceU nucleic acid in which the 
Intent is no, ordinarily found. Exogenous DNA segments are expressed to yteld 
exogenous poises. ^ ^ ^ , ^ ^ ^ ^ rf dn a ^ 

■ ,« , „ Thus eenes include coding sequences and/or the regulatory 
with a hiologtca, functton Thus, genes in „ onexpre ssed DNA segments 

sequences required for their expression. Genes also 
that for example, form recognition sequences for outer protems. 

The term "wild-type" means that the nucleic acid fragment does no, 
comprise any nations. A "wild-type" protein means tha, tire pro,=in will be active at a 
ZZ Z ity found in nature and typically will comprise the amino acid sequence 
I! 1 naturl ,n an aspect, the term "wiid ,pe" or "parent sequence" can indicate a 
startingorreferencesequencepnortoamanipulationoftheinventton. 

"Substantially pure" means an object speces is the predomman, speces 
» oresen, (i e on a molar basis i, is more abundantman any outer indivtdual 
3 present u.c. on „„f„ a blv a substantially purified 

macromolecular species in the composition), and preferably 

fraction is a composition wherein the object species compnses a, leas, abou, 50 p rcent 
I a, basis) of al, macromolecular speces present. Generally, a substantia y pure 
opposition wi« comprise more than abou. SO to 90 percent of aU -romolecul 
« spec es present in the composition. Most preferably, the object speces ,s purified , 

— nal ^^-^^^1^ 

macromolecular species. Solvent species, small molecules ( 

ion species are not considered macromolecular species. 
.„ Percentage sequence identity is calculated by comparing two optimally 

aligned sequences over "the window of comparison, determining tine number of positions 
wlh t e identical nucleic add base occurs ,n both sequences ,o yiel ,he number of 
mlched positions, dividing the number of matched positions by the tota! number of 



PCT/US98/25698 



WO 99/29902 



10 



n«timal alignment of sequences for aligning a 
Reiease 7 0 Genetics Computer Group, 575 Sctence Dr., Madtson. Wl. 

po.ypepude or ^nucleotide — ^d « I L no, L ta— -By modifted 
can be tsorated from a source ,n nature »d w*ch h ^ ^ 

rCattonsnip witb ^ = - ^ ~ * £ ^ scrip l of ths codmg 

op erab, y UnKed to a codi» g ^ " ^ bring are typicaUy 

::r;re^::::,— — 

promoter by several abases and inuonic sequences may be of vanaHe , 

, onerablv linked but not contiguous. 

in the human genome, the human CD4 gene g 

„o «f these two genes indicate that they are nigmy 
?S since the sequences and structures of these two gen 

boltoj and bo.b g enes encode a pro.ern whicb ftmcons ,n , g naHn g T ce„ 

• u u uur class Il-restricted antigen recognition, 
activate tnrou^b MHC c£^_ ^ ^ ^ ^ a ^ by ^ ^ 

• strands derived from the different parental duplex 
between complementary ^■""*£„ .0 double-stranded DMA generated bv 
,0 modules, whereas the ,enn homodup e> e - ^ 
base pairing between complementary s,n g le strands 



duplex molecules. 
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Tta term duple* ™ A » "« 3 

"'° OP " ZZr vartan, seance is a se q uence snowm, suhstanUal vanat,on 

sequence at one or more positions. 

DmJ IFP n^r.RTPTIQN 

I. General _ . polvnuc leotide toward 

The invention provides methods ot evolving v . 

,eas, two polynucleotide variant sequences — 1J S a. annealed 

in vuro at the regtons of smulanty. ^ ^ 

form heteroduplexes, in wh,ch component str^ds o g ^ 

products of annealing ate exposed to enzymes of a DNA rep P 
Ipiicatton system, tltat repatrs — d P-^^ QNA repair 

° — — ptod^ateexposed 
svstem Alternatively, exposure can be m vitro. 

system, a ,i n nnl DNA repair systems. Exposure of 

to cellular extracts containing functional DNA repa y 

— exestoaDNArep^^ 

heterodupiexes due ,0 DNA m.smatchmg. Th repan T ^ 
25 recomW „an on.npromot.ng, — ' 
DNA tepa.t ptocess ,s typ.c lly effected ^ ^ ^ ^ 

mol ecu,e and a, any part.cular ^ polymlcUotides 

paired. The testing popuiauon can thus conta 
encompassing an essentiaUy random teassorunent o pom s o 

. The rtooulation of recombinant polynucleotides is then sere 
30 parental strands. The popuiauon aorop er<v of the polvnucleotide 

expression product thereof, such as mRNA or a protein. 
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n SuhstTatcs Fnr Shuffling t -A„.y a . 

Tfc c substra.es for shuffling are variants of a reference polynucieottde that 
sh ow some region(s) of similanty wtth the reference and Cher reg,on(s> or potnt.s, of 
d,verg=nce. Reg.ons of similariry should be suffrciem .o support anncalmg of 
po.ynucleotides such that s.ab>e heteroduplexes can be formed. Variants 

1 —a, seance identity with each other (e.g., at >eas, 50%, 75%. *>% or 
mere should be a, least suffice diversity between substrates tha, recombmatton .can 
g e„era,e more diverse product, than mere are starting materia,. Thus, mere must be , 
L, two substrates dtffering in a, teas, two positions. The degree of dtverstty epends on 
the length of the substrate being recombined and the extent of the funcuonai change to be 
evoived. Diversity a, between 0,-25% of positions ,s typ.cal. Recombmauon of 
mutations from very cioseiy related genes or even whole sections of seances from more 
distandy related genes or sets of genes can enhance the rate of evolutton and the 

cquisi on of desirable new properties. Recombinant to create chunertc or mosa. 
ZL be usefu, in order ,0 combine desirable features of two or more parents mto a 
Z g -e or set of genes, or to create novel feamres not found m the ^parents. 

2 number of different substrates ,0 be combined can vary widely in « from two to 
10 100 1000. to more than 10\ 10'- or 10' members. 

The initial small population of the specifc nucleic acid sconces havtng 
nations may be created by a number o, different method, Mutations may be c a,d 
by error-prone PCR. Error-prone PCR uses .ow-fidelity polymertzatton condt.tons to 
Loduce a low level of point — randomly over a long seance. Al,crnauv.y 
Itations can be mtroduced into the template polynucleotide by oligonucleo „de-d,rected 
5 mutagens, In oligonucleotide-directed mutagenesis, a short sequence of the 

polynucleotrde is removed from the polynucleotide using restncon enzyme d.ges ton^ 
and is replaced with a synthettc poiynucleotide in which vanous bases have been altered 
:l t e ortgina, sequence. The P o,ynuc,eo,de sequence can also be ahered by chemtca, 
L s enes,s Chemtca, mutagens tnclude. for example, sod.um bisu, ,e. nttrous actd. 
30 hydroxviamme. hydrazine or form.c acd. Other agents wh,ch are analogues of 
Leoude precursors .nc.ude nitrosoguamd.ne. 5-bromouracil. 2 am.nopunne. or 
acridme. GeneraHy. these agents are added to the PCR reacon ,n place of the nucleoude 
precursor therebv mutattng the sequence. Intercalate agents such as proflav.ne. 
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acrifiavine. quinacnne and the like can also be used. Random mutageneses of the 
polynucleotide sequence can also be achieved by Nation with X-rays or ultraviolet 
light. Generally, plasrmd DNA or DNA fragments so mutagenized are introduced mto E. 
coli and propagated as a pool or library of mutant plasmids. 

Alternatively the small mixed population of specific nucleic acids can be 
found in nature in the form of different alleles of the same gene or the same gene from 
different related species (i.e., cognate genes). Alternatively, substrates can be related but 
nonallelic genes, such as the immunoglobulin genes. Diversity can also be the result of 
previous recombination or shuffling. Diversity can also result from resynthesizmg genes 
encoding natural proteins with alternative codon usage. 

The starting substrates encode variant forms of sequences to be evolved. 
In some methods, the substrates encode variant forms of a protein for which evolution of 
a new or modified property is desired. In other methods, the substrates can encode 
variant forms of a plurality of genes constituting a multigene pathway. In such methods, 
variation can occur in one or any number of the component genes. In other methods, 
substrates can contain variants segments to be evolved as DNA or RN A binding 
sequences. In methods, in which starting substrates containing coding sequences, any 
essential regulatory sequences, such as a promoter and polyadenylation sequence, 
required for expression may also be present as a component of the substrate. 
Alternatively, such regulatory sequences can be provided as components of vectors used 

for cloning the substrates. 

The starting substrates can vary in length from about 50, 250, 1000, 
10 000 100,000, 10 6 or more bases. The starting substrates can be provided in double- or 
single-stranded form. The starting substrates can be DNA or RNA and analogs thereof. 
If DNA, the starting substrates can be genomic or cDNA. If the substrates are RNA. the 
substrates are typically reverse-transcribed to cDNA before heteroduplex formation. 
Substrates can be provided as cloned fragments, chemically synthesized fragments or 
PGR amplification products. Substrates can derive from chromosomal, plasrmd or viral 
sources. In some methods, substrates are provided in concatemeric form. 



30 



TIT Procedure for fierier -""" H^eroduplexes 

Heteroduplexes are generated from double stranded DNA substrates, by 
denaturing the DNA substrates and incubating under annealing conditions. Hybridization 



W O 99/29902 



PCT/US98/25698 



12- 



condi.ions fo, heteroduplex tendon are ********** and « diff -" » 

circumstances. Longer sequences hybridize specifically a. higher .empera^re. 
Generally, hybridiza,on condLions m selecud to be about 25'C lower than <he .hermal 
melon, poin, (Tm) for .he specific sequence a. a defined .cue suenph and pH. The Tm 
; , s the .ensure (under defined ionic strength. pH. and nucleic acd concentration) a, 
which 50% of .he probes complement ,0 ,he .arge. sequence hybridize .o .he .arge. 

sequence at equilibrium. 

Exemplary conditions for dena.ura.ion and renaturation of double stranded 

substrates are as foUows. Equimolar concentrations (- !.0 - 5.0 nM, ' » 

0 mixed in 1 x SSPE buffer (IS. mM NaCl. .0 mM EDTA, 10 mMNarWO.. pH 7.4, 
After heating a, 96°C for 1 0 minutes, the reaction mixture is immedtately cooled a. 
for 5 minutes, The m,x.ure is then incubated a, 68'C for 2-6 nr. Dena.ura.ton and 
sealing can also be carried „u, by ,he addhion and remova, of a denaturant such as 
NaOH. The process is me same for single stranded DNA substrates, exeep. that the 
,5 denaturing step may be omitted for short sequences. 

By appropriate design of substrates for heteroduplex formauon. ,t rs 
possible to achieve selection for heteroduplexes relative .0 reformed parental 
homoduplexes. Homoduplexes merely reconstruct parental sub*ra,es and effectively 
dilute recombinant products in subsequen, screening s.e P s. In general, seieetton ,s 
,„ achieved bv designing subs.ra.es such ma. he,=r„dup.exes are formed in open-crcles, 

whereas homoduplexes are formed as linear molecules. A subsequen. transformation step 
results in substantial enrichmen, (e.g., 100-fold) to, ,he circular heteroduplexes. 

Figure 1 shows a method in which two substrate sequences in separate 
vectors are PCR-amplified using .wo different sets of primers (PI, P2 and P3. P4). 
,5 Typically, firs, and second substrates are inserted into separate copies of the same vector. 
" The .wo differen, pairs of pnmers initiate amplification a. different po.n.s on ,he .wo 
vecors. Fig. . shows an arrangement in which the P1/P2 primer pairs initiates 
amplication at one of the two boundaries of the vector with .he substrate and ,he P„P2 
pr ,mer pair initiates replication a, the other boundary ,n a second vector. The two prtmers 
„ each pnmcr pair prime amplifica.ion in opposi.e directions around a circular plasmid. 

, a u„ amnlification are double-stranded linearized 
The amplification products generated by this amplication 

ector molecules in which the first and second substrates occur at opposite ends of the 
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vector. The amplification products are mixed, denatured and annealed. Mixing and 
denaturation can be performed in either order. Reannealing generates two linear 
homoduplexes. and an open circular heteroduplex containing one nick in each strand, at 
the initiation point of PCR amplification. Introduction of the amplification products into 
5 host cells selects for the heteroduplexes relative to the homoduplexes because the former 
transform much more efficiently than the latter. 

It is not essential in the above scheme that amplification is initiated at the 
interface between substrate and the rest of the vector. Rather, amplification can be 
initiated at any points on two vectors bearing substrates provided that the amplification is 
1 0 initiated at different points between the vectors. In the general case, such amplification 
generates two linearized vectors in which the first and second substrates respectively 
occupy different positions relative to the remainder of the vector. Denaturation and 
reannealing generator heteroduplexes similar to that shown in Fig. 1 , except that the nicks 
occur within the vector component rather than at the interface between plasmid and 
1 5 substrate. Initiation of amplification outside the substrate component of a vector has the 
advantage that it is not necessary to design primers specific for the substrate borne by the 
vector. 

Although Fig. 1 is exemplified for two substrates, the above scheme can be 
extended to any number of substrates. For example, an initial population of vector 
20 bearing substrates can be divided into two pools. One pool is PCR-amplified from one 

set of primers, and the other pool from another. The amplification products are denatured 
and annealed as before. Heteroduplexes can form containing one strand from any 
substrate in the first pool and one strand from any substrate in the second pool. 
Alternatively, three or more substrates cloned into multiple copies of a vector can be 
25 subjected to amplification with amplification in each vector starting at a different point. 
For each substrate, this process generates amplification products varying in how flanking 
vector DNA is divided on the two sides of the substrate. For example, one amplification 
product might have most of the vector on one side of the substrate, another amplification 
product might have most of the vector on the other side of the substrate, and a further 
30 amplification product might have an equal division of vector sequence flanking the 
substrate. In the subsequent annealing step, a strand of substrate can form a circular 
heteroduplex with a strand of any other substrate, but strands of the same substrate can 
only reanneal with each other to form a linear homoduplex. In a still further variation. 
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multiple substrates can be performed by performing multiple iterations of the scheme in 
Fig. 1 . After the first iteration, recombinant polynucleotides in a vector, undergo 
heteroduplex formation with a third substrate incorporated into a further copy of the 
vector. The vector bearing the recombinant polynucleotides and the vector bearing the 
third substrate are separately PCR amplified from different primer pairs. The 
amplification products are then denatured and annealed. The process can be repeated 
further times to allow recombination with further substrates. 

An alternative scheme for heteroduplex formation is shown in Fig. 2. 
Here, first and second substrates are incorporated into separate copies of a vector. The 
two copies are then respectively digested with different restriction enzymes. Fig. 2 shows 
an arrangement in which, the restriction enzymes cut at opposite boundaries between 
substrates and vector, but all that is necessary is.to use two different restriction enzymes 
that cut at different places. Digestion generates linearized first and second vector bearing 
first and second substrates, the first and second substrates occupying different positions 
relative to the remaining vector sequences. Denaturation and reannealing generates open 
circular heteroduplexes and linear homoduplexes. The scheme can be extended to 
recombination between more than two substrates using analogous strategies to those 
described with respect to Fig. 1. In one variation, two pools of substrates are formed, and 
each is separately cloned into vector. The two pools are then cute with different enzymes, 
and annealing proceeds as for two substrates. In another variation, three or more 
substrates can be cloned into three or more copies of vector, and the three or more result 
molecules cut with three or more enzymes, cutting at three or more sites. This generates 
three different linearized vector forms differing in the division of vector sequences 
flanking the substrate moiety in the vectors. Alternatively, any number of substrates can 
be recombined pairwise in an iterative fashion with products of one round of 
recombination annealing with a fresh substrate in each round. 

In a further variation, heteroduplexes can be formed from substrates 
molecules in vector-free form, and the heteroduplexes subsequently cloned into vectors. 
Such can be achieved by asymmetric amplification of first and second substrates as 
shown in Fig. 3. Asymmetric or single primer PCR amplifies only one strand of a duplex. 
By appropriate selection of primers, opposite strands can be amplified from two different 
substrates. On reannealing amplification products, heteroduplexes are formed from 
opposite strands of the two substrates. Because only one strand is amplified from each 
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substrate, reannealing does not reform homoduplexes (other than for small quantities of 
unamplified substrate). The process can be extended to allow recombination of any 
number of substrates using analogous strategies to those described with respect to Fig. 1. 
For example, substrates can be divided into two pools, and each pool subject to the same 
5 asymmetric amplification, such that amplification products of one pool can only anneal 
with amplification products of the other pool , and not with each other. Alternatively, 
shuffling can proceed pairwise in an iterative manner, in which recombinants formed 
from heteroduplexes of first and second substrates, are subsequently subjected to 
heteroduplex formation with a third substrate. Point mutations can also be introduced at a 
10 desired level during PCR amplification. 

Fig. 4 shows another approach of selecting for heteroduplexes relative to 
homoduplexes. First and second substrates are isolated by PCR amplification from 
separate vectors. The substrates are denatured and allowed to anneal forming both 
heteroduplexes and reconstructed homoduplexes. The products of annealing are digested 
1 5 with restriction enzymes X and Y. X has a site in the first substrate but not the second 
substrate, and vice versa for Y. Enzyme X cuts reconstructed homoduplex from the first 
substrate and enzyme Y cuts reconstructed homoduplex from the second substrate. 
Neither enzyme cuts heteroduplexes. Heteroduplexes can effectively be separated from 
restriction fragments of homoduplexes by further cleavage with enzymes A and B having 
20 sites proximate to the ends of both the first and second substrates, and ligation of the 

products into vector having cohesive ends compatible with ends resulting from digestion 
with A and B. Only heteroduplexes cut with A and B can ligate with the vector. 
Alternatively, heteroduplexes can be separated from restriction fragments of 
homoduplexes by size selection on gels. The above process can be generalized to N 
25 substrates by cleaving the mixture of heteroduplexes and homoduplexes with N enzymes, 
each one of which cuts a different substrate and no other substrate. Heteroduplexes can 
be formed by directional cloning. Two substrates for heteroduplex formation can be 
obtained by PCR amplification of chromosomal DNA and joined to opposite ends of a 
linear vector. Directional cloning can be achieved by digesting the vector with two 
30 different enzymes, and digesting or adapting first and second substrates to be respectively 
compatible with cohesive ends of only of the two enzymes used to cut the vector. The 
first and second substrates can thus be ligated at opposite ends of a linearized vector 
fragment. This scheme can be extended to any number of substrates by using principles 
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analogous to those described for Fig. 1 . For example, substrates can be divided into two 
pools before ligation to the vector. Alternatively, recombinant products formed by 
heteroduplex formation of first and second substrates, can subsequently undergo 
heteroduplex formation with a third substrate. 

5 

IV. Vectors and Transformation 

In general, substrates are incorporated into vectors either before or after 
the heteroduplex formation step. A variety of cloning vectors typically used in genetic 
engineering are suitable. 
0 The vectors containing the DNA segments of interest can be transferred 

into the host cell by standard methods, depending on the type of cellular host. For 
example, calcium chloride transformation is commonly utilized for prokaryotic cells, 
whereas calcium phosphate treatment. Lipofection, or electroporation may be used for 
other cellular hosts. Other methods used to transform mammalian cells include the use of 
5 Polybrene, protoplast fusion, liposomes, electroporation, and microinjection, and 

biolisitics (see, generally, Sambrook et al., supra). Viral vectors can also be packaged in 
vitro and introduced by infection. The choice of vector depends on the host cells. In 
general, a suitable vector has an origin of replication recognized in the desired host cell, a 
selection maker capable of being expressed in the intended host cells and/or regulatory 
20 sequences to support expression of genes within substrates being shuffled. 

V Types of Host Cells 

In general any type of cells supporting DNA repair and replication of 
heteroduplexes introduced into the cells can be used. Cells of particular interest are the 

25 standard cell types commonly used in genetic engineering, such as bacteria, particularly, 
E. coli (16, 17). Suitable E. coli strains include E. coli mutS, mutU dam , and/or recA* . 
E.coli XL-10-Gold ([Tet r A(mcrA)183 A(mcrCB~hsdSMR-mrr) 1 73 endAl supE44 thi-1 
recAl gyrA96 relAl lac Hte] [F' proAB lacPZAMIS TnlO CTef) Amy Cam J), E. coli 
ES1301 mutS [Genotype: lacZ53. mu(S201 : :Tn5, thyA36. rha-5, meiBL deoC. INfrrnD- 

30 rrnE)] (20. 24. 28-42). Preferred E. coli strains are E.coli SCSI 10 [Genotype: rpsi 
(Sir), thr. leu. enda. thi-1. lacy. galk. gait, ara tona. tsx, dam. dcm. supE44. A(lac- 
proAB), [F. traD36. proA~ B~ lacVZAMl 5\ which have normal cellular mismatch repair 
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systems (17). This strain type repairs mismatches and unmatches in the heteroduplex 
with little strand-specific preference. Further, because this strain is darn and dcm\ 
piasmid isolated from the strain is unmethylated and therefore particularly amenable for 
further rounds of DNA duplex formation/mismatch repair (see below). Other suitable 
5 bacterial cells include gram-negative and gram-positive, such as Bacillus. Pseudomonas, 
and Salmonella. 

Eukaryotic organisms are also able to carry out mismatch repair (43-48). 
Mismatch repair systems in both prokaryotes and eukaryotes are thought to play an 
important role in the maintenance of genetic fidelity during DNA replication. Some of 
10 the genes that play important roles in mismatch repair in prokaryotes, particularly mutS 

and mutL. have homologs in eukaryotes. in the outcome of genetic recombinations, and in 
genome stability. Wild-type or mutant S. cerevisiae has been shown to carry out 
mismatch repair of heteroduplexes (49-56), as have COS-1 monkey cells (57). Preferred 
strains of yeast are Picchia and Saccharomyces. Mammalian cells have been shown to 
1 5 have the capacity to repair G-T to G-C base pairs by a short-patch mechanism (38,58- 

63). Mammalian cells (e.g., mouse, hamster, primate, human), both cell lines and primary 
cultures can also be used. Such cells include stem cells, including embryonic stem cells, 
zygotes, fibroblasts, lymphocytes, Chinese hamster ovary (CHO), mouse fibroblasts 
(NIH3T3), kidney, liver, muscle, and skin cells. Other eucaryotic cells of interest include 
20 plant cells, such as maize, rice, wheat, cotton, soybean, sugarcane, tobacco, and 

arabidopsis; fish, algae, fungi (aspergillus, podospora. neurospora), insect (e.g.. baculo 
lepidoptera) (see, Winnacker. "From Genes to Clones. ,, VCH Publishers. N.Y., (1987), 
which is incorporated herein by reference). 

In vivo repair occurs in a wide variety of prokaryotic and eukaryotic cells. 
25 Use of mammalian cells is advantage in certain application in which substrates encode 

polypeptides that are expressed only in mammalian cells or which are intended for use in 
mammalian cells. However, bacterial and yeast cells are advantageous for screening 
large libraries due to the higher transformation frequencies attainable in these strains. 

30 V. In Vitro DNA Repair Systems 

As an alternative to introducing annealed products into host cells, annealed 
products can be exposed a DNA repair system in vitro. The DNA repair system can be 
obtained as extracts from repair-competent E. coli. yeast or any other cells (64-67). 
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Repair-competent cells are lysed in appropnate buffer and supplemented with 

Repair comp lrans formed into competent cells 

nucleotides. DNA is incubated in this cell extraci 
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bind a given protein target, and capacity to regulate expression of operably linked coding 
sequences. Some of the above properties, such as drug resistance, can be selected by 
plating cells on the drug. Other properties, such as the influence of a regulatory sequence 
on expression, can be screened by detecting appearance of the expression product of a 
5 reporter gene linked to the regulatory sequence. Other properties, such as capacity of an 
expressed protein to be secreted, can be screened by FACS™, using a labelled antibody to 
the protein. Other properties, such as immunogenicity or lack thereof, can be screened by 
isolating protein from individual cells or pools of cells, and analyzing the protein in vitro 
or in a laboratory animal. 

10 

VII. Variations 

1. Demethvlation 

Most cell types methylate DNA in some manner, with the pattern of 
methylation differing between cells types. Sites of methylation include 5-methylcytosine 
1 5 (m 5 C), N4-methylcytosine (m 4 C) and N 6 -methyladenine (m 6 A), 5 -hydroxy me thy Icytosine 
(hm 5 C) and 5-hydroxymethyluracil (hm 5 U). In E. coli, methylation is effected by Dam 
and Dcm enzymes. The methylase specified by the dam gene methylates the N6-position 
of the adenine residue in the sequence GATC, and the methylase specified by the dcm 
gene methylates the C5-position of the internal cytosine residue in the sequence 
20 CCWGG. DNA from plants and mammal is often subject to CG methylation meaning 
that CG or CNG sequences are methylated. Possible effects of methylated on cellular 
repair are discussed by references 1 8-20. 

In some methods, DNA substrates for heteroduplex formation are at least 
partially demethylated on one or both strands, preferably the latter. Demethyiation of 
25 substrate DNA promotes efficient and random repair of the heteroduplexes. In 

heteroduplexes formed with one strand dam-methylated and one strand unmethylated, 
repair is biased to the unmethylated strand, with the methylated strand serving as the 
template for correction. If neither strand is methylated, mismatch repair occurrs. but 
showes insignificant strand preference (23. 24). 
30 Demethyiation can be performed in a variety of ways. In some methods, 

substrate DNA is demethylated by PCR-amplification. In some instances. DNA 
demethyiation is accomplished in one of the PCR steps in the heteroduplex formation 
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procedures described above. In other methods, an additional PCR step is performed to 
effect demethylation. In other methods, demethylation is effected by passaging substrate 
DNA through methylation deficient host cells (e.g. an E. coli dam'dcm strain). In other 
methods, substrate DNA is demethylated in vitro using a demethylating enzymes. 
Demethylated DNA is used for heteroduplex formation using the same procedures 
described above. Heteroduplexes are subsequently introduced into DNA-repair-proficient 
but restriction-enzyme-defective cells to prevent degradation of the unmethylated 
heteroduplexes. 

2. Sealing Nicks 

Several of the methods for heteroduplex formation described above result 
in circular heteroduplexes bearing nicks in each strand. These nicks can be sealed before 
introducing heteroduplexes into host cells. Sealing can be effected by treatment with 
DNA ligase under standard ligating conditions. Ligation forms a phosphodiester bond to 
link two adjacent bases separated by a nick in one strand of double helix of DNA. 
Sealing of nicks increases the frequency of recombination after introduction of 
heteroduplexes into host cells. 

3. Error Prone PCR Attendant To Amplification 

Several of the formats described above include a PCR amplification step. 
Optionally, such a step can be performed under mutagenic conditions to induce additional 
diversity between substrates. 

VIII. Other Shuffling Methods 
5 The methods of heteroduplex formation described above can be used in 

conjunction with other shuffling methods. For example, one can perform one cycle of 
heteroduplex shuffling, screening or selection, followed by a cycle of shuffling by another 
method, followed by a further cycle of screening or selection. Other shuffling formats are 
described by WO 95/22625: US 5,605.793; US 5.81 1.238; WO 96/19256; Stemmer. 
0 Science 270. 1510 (1995); Stemmer et aL Gene. 164.49-53 (1995); Stemmer. 

Bio/Technology, 13. 549-553 (1995); Stemmer. Proc. Natl. Acad Sci. USA 91. 10747- 
10751 (1994); Stemmer. Nature 370. 389-391 (1994); Crameri et al.. Nature Medicine. 
2(1): 1-3. (1996): Crameri et al.. Nature Biotechnology 14. 315-319 (1996); WO 98/42727; 
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WO 98/41622; WO 98/05764 and WO 98/42728. WO 98/27230 (each of which is 
incorporated by reference in its entirety for all purposes). 



IX. Protein Analogs 

5 Proteins isolated by the methods also serve as lead compounds for the 

development of derivative compounds. The derivative compounds can include chemical 
modifications of amino acids or replace amino acids with chemical structures. The 
analogs should have a stabilized electronic configuration and molecular conformation that 
allows key functional groups to be presented in substantially the same way as a lead 
10 protein. In particular, the non-peptic compounds have spatial electronic properties which 
are comparable to the polypeptide binding region, but will typically be much smaller 
molecules than the polypeptides, frequently having a molecular weight below about 2 
CHD and preferably below about 1 CHD. Identification of such non-peptic compounds 
can be performed through several standard methods such as self-consistent field (CSF) 
15 analysis, configuration interaction (CHI) analysis, and normal mode dynamics analysis. 
Computer programs for implementing these techniques are readily available. See Rein et 
al., Computer-Assisted Modeling of Receptor-Ligand Interactions (Alan Liss, New York, 
1989). 

20 IX. Pharmaceutical Compositions 

Polynucleotides, their expression products, and secondary metabolites 
whose formation is catalyzed by expression products, generated by the above methods are 
optionally formulated as pharmaceutical compositions. Such a composition comprises 
one or more active agents, and a pharmaceutically acceptable carrier. A variety of 

25 aqueous carriers can be used, e.g., water, buffered water, phosphate-buffered saline 

(PBS), 0.4% saline, 0.3% glycine, human albumin solution and the like. These solutions 
are sterile and generally free of particulate matter. The compositions may contain 
pharmaceutically acceptable auxiliary substances as required to approximate physio- 
logical conditions such as pH adjusting and buffering agents, toxicity adjusting agents 

30 and the like, for example, sodium acetate, sodium chloride, potassium chloride, calcium 
chloride and sodium is selected primarily based on fluid volumes, viscosities, and so 
forth, in accordance with the particular mode of administration selected. 
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EXAMPLES 

EXAMPLE 1 . Novel Rhizobium Flag Genes From Recombination Of Rhizobium Lupini 
Flag And Rhizobium Meliloti FlaA 

Bacterial flagella have a helical filament, a proximal hook and a basal 
5 body with the flagellar motor (68). This basic design has been extensively examined in 
£. coli and S. typhimurium and is broadly applicable to many other bacteria as well as 
some archaea. The long helical filaments are polymers assembled from flagellin subunits, 
whose molecular weights range between 20,000 and 65,000, depending on the bacterial 
species (69). Two types of flagellar filaments, named plain and complex, have been 
10 distinguished by their electron microscopically determined surface structures (70). Plain 
filaments have a smooth surface with faint helical lines, whereas complex filaments 
exhibit a conspicuous helical pattern of alternating ridges and grooves. These 
characteristics of complex flagellar filaments are considered to be responsible for the 
brittle and (by implication) rigid structure that enables them to propel bacteria efficiently 
1 5 in viscous media (71-73). Whereas flagella with plain filaments can alternate between 

clockwise and counter clockwise rotation (68), all known flagella with complex filaments 
rotate only clockwise with intermittent stops (74). Since this latter navigation pattern is 
found throughout bacteria and archaea, it has been suggested that complex flagella may 
reflect the common background of an ancient, basic motility design (69). 
20 Differing from plain bacterial flagella in the fine structure of their 

filaments dominated by conspicuous helical bands and in their fragility, the filaments are 
also resistant against heat decomposition (72). Schmitt et al. (75) showed that 
bacteriophage 7-7-1 specifically adsorbs to the complex flagella of RAupini H 13-3 and 
requires motility for a productive infection of its host. Though the flagellins from R. 
25 meliloti and R. lupini are quite similar, bacteriophage 7-7-1 does not infect R. meliloti. 
Until now complex flagella have been observed in only three species of soil bacteria: 
Pseudomonas rhodos (73), R.meliloti (76), and RAupini HI 3-3 (70, 72). Cells of RAupini 
HI 3-3 posses 5 to 10 peritrichously inserted complex flagella. which were first isolated 
and analyzed by high resolution electron microscopy and by optical diffraction (70). 
30 Maruyama et al. (77) further found that a higner content of hydrophobic 

amino acid residues in the complex filament may be one of the main reasons for the 
unusual properties of complex flagella. By measuring mass per unit length and obtaining 
three-dimensional reconstruction from electron micrographs. Trachtenberc et al. (73. 78) 
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suggested that the complex filaments of R. lupini are composed of functional dimers. 
Figure 6 shows the comparison between the deduced amino acid sequence of the R. lupini 
HI 3-3 FlaA and the deduced amino acid sequence of the R. meliloti FlaA. Perfect 
matches are indicated by vertical lines, and conservative exchanges are indicated by 
5 colons. The overall identity is 56%. The RAupiniflaA and R.meliloti JlaA were subjected 
to in vitro heteroduplex formation followed by in vivo repair in order to create novel 
FlaA molecules and structures. 

A. Methods 

1 0 pRL20 containing R. lupini HI 3-3 JlaA gene and pRM40 containing 

R.meliloti JlaA gene are shown in Figs. 6 A and 6B. These plasmids were isolated from £ 
coli SCSI 10 (free from dam- and dcm-type methylation). 

About 3.0 pg. of unmethylated pRL20 and pRM40 DNA were digested with Bam HI and 
Eco RL respectively, at 37°C for 1 hour. After agarose gel separation, the linearized 

1 5 DNA was purified with Wizard PCR Prep kit (Promega, WI, USA). 

Equimolar concentrations (2.5 nM) of the linearized unmethylated pRL20 and pRM40 
were mixed in 1 x SSPE buffer (180 mMNaCl, 1 mM EDTA, 10 mM NaH2P04, pH 
7.4). After heating at 96°C for 10 minutes, the reaction mixture was immediately cooled 
at 0°C for 5 minutes. The mixture was incubated at 68°C for 2 hour for heteroduplexes to 

20 form. 

One microliter of the reaction mixture was used to transform 50 \x\ of E. 
coli ES 1301 mutS, £ coli SCSI 10 and £ coli JM109 competent cells. The 
transformation efficiency with E. coli JM109 competent cells was about seven times 
higher than that of E. coli SCSI 10 and ten times higher than that of £ coli ES1301 mutS, 
25 although the overall transformation efficiencies were 1 0-200 times lower than those of 
control transformations with the close, covalent and circular pUC19 plasmid. 

Two clones were selected at random from the E. coli SCSI 10 
transformants and two from £ coli ESI 301 mutS transformants. and plasmid DNA was 
isolated from these four clones for further DNA sequencing analysis. 
30 B. Results 

Figure 7 shows (a) the sequence of SCS01 (clone#l from £ coli SCSI 10 
transformant library;, (b) the sequence of SCS02 (clone #2 from £ coli SCSI 10 
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transformant library), (c) the sequence of ES01 (clone #1 from E. coli ES1301 
transformant library), and (d) the sequence of ES02 (clone #2 from E. coli ESI 301 
transformant library). All four sequences were different from wild-type R. lupiniflaA and 
R. melilotiflaA sequences. Clones SCS02. ES01 and ES02 all contain a complete open- 
5 reading frame, but SCS01 was truncated. Figure 8 shows that recombination mainly 
occurred in the loop regions (unmatched regions). The flaA mutant library generated 
from R. melilotiflaA and R. lupiniflaA can be transformed into E. coli SCSI 10. ESI 301, 
XLIO-Gold and JM109, and transformants screened for functional FlaA recombinants. 

10 EXAMPLE 2. Directed Evolution Of ECB Deacvlase For Variants With Enhanced 
Specific Activity 

Streptomyces are among the most important industrial microorganisms due 
to their ability to produce numerous important secondary metabolites (including many 
antibiotics) as well as large amounts of enzymes. The approach described here can be 
15 used with little modification for directed evolution of native Streptomyces enzymes, 

some or all of the genes in a metabolic pathways, as well as other heterologous enzymes 
expressed in Streptomyces, 

New antifungal agents are critically needed by the large and growing 
numbers of immune-compromised AIDS, organ transplant and cancer chemotherapy 
20 patients who suffer opportunistic infections. Echinocandin B (ECB), a lipopeptide 
produced by some species of Aspergillus, has been studied extensively as a potential 
antifungal. Various antifungal agents with significantly reduced toxicity have been 
generated by replacing the linoleic acid side chain of A. nidulans echinocandin B with 
different aryl side chains (79-83). The cyclic hexapeptide ECB nucleus precursor for the 
25 chemical acylation is obtained by enzymatic hydrolysis of ECB using Actinoplanes 

utahensis ECB deacvlase. To maximize the conversion of ECB into intact nucleus, this 
reaction is carried out at pH 5.5 with a small amount of miscible organic solvent to 
solubilize the ECB substrate. The product cyclic hexapeptide nucleus is unstable at pH 
above 5.5 during the long incubation required to fully deacylate ECB (84). The pH 
30 optimum of ECB deacvlase. however, is 8.0-8.5 and its activity is reduced at pH 5.5 and 
in the presence of more than 2.5% ethanol (84). To improve production of ECB nucleus 
it is necessary to increase the activity of the ECB deacvlase under these process-relevant 
conditions. 
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Relatively little is known about ECB deacylase. The enzyme is a 
heterodimer whose two subunits are derived by processing of a single precursor protein 
(83). The 19.9 kD a-subunit is separated from the 60.4 kD p-subunit by a 15-amino acid 
spacer peptide that is removed along with a signal peptide and another spacer peptide in 
5 the native organism. The polypeptide is also expressed and processed into functional 

enzyme in Streptomyces lividans* the organism used for large-scale conversion of ECB by 
recombinant ECB deacylase. The three-dimensional structure of the enzyme has not been 
determined, and its sequence shows so little similarity to other possibly related enzymes 
such as penicillin acylase that a structural model reliable enough to guide a rational effort 
1 0 to engineer the ECB deacylase will be difficult to build. We therefore decided to use 
directed evolution (85) to improve this important activity. 

Protocols suitable for mutagenic PCR and random-priming recombination 
of the 2.4 kb ECB deacylase gene (73% G+C) have been described recently (86). Here, 
we further describe the use of heteroduplex recombination to generate new ECB 
15 deacylase with enhanced specific activity. 

In this case, two Actinoplanes ulahensis ECB deacylase mutants, M7-2 
and Ml 6, which show higher specific activity at pH 5.5 and in the presence of 10% 
MeOH were recombined using technique of the in vitro heteroduplex formation and in 
vivo mismatch repair . 

20 Figure 12 shows the physical maps of plasmids pM7-2 and pM16 which 

contain the genes for the M7-2 and Ml 6 ECB deacylase mutants. Mutant M7-2 was 
obtained through mutagenic PCR performed directly on whole Streptomyces lividans cells 
containing wild-type ECB deacylase gene, expressed from plasmid pSHP 150-2*. 
Streptomyces with pM7-2 show 1.5 times the specific activity of cells expressing the 

25 wild-type ECB deacylase (86). Clone pM16 was obtained using the random-priming 

recombination technique as described (86, 87). It shows 2.4 times specific activity of the 
wild-type ECB deacylase clone. 

A. Methods : 

30 M7-2 and Ml 6 plasmid DNA (pM7-2 and pM16) (Fig. 9) were purified 

from £. coli SCSI 10 (in separate reactions). About 5.0 \x% of unmethylated M7-2 and 
Ml 6 DNA were digested with XJio I and Psh AI. respectively, at 37°C for 1 hour (Fig. 
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10). After agarose gel separation, the linearized DNA was purified using a Wizard PCR 
Prep Kit (Promega, WL USA). 

Equimolar concentrations (2.0 nM) of the linearized unmethylated pM7-2 and pM16 
DNA were mixed in 1 x SSPE buffer (lx SSPE: 180 mM NaCl, 1.0 mM EDTA, 10 mM 
5 NaH 2 P0 4 , pH 7.4). After heating at 96 °C for 10 minutes, the reaction mixture is 
immediately cooled at 0 °C for 5 minutes. The mixture was incubated at 68 °C for 3 
hours to promote formation of heteroduplexes. 

One microliter of the reaction mixture was used to transform 50 ^1 of 
E.coli ES1301 mutS, SCSI 10 and JM109 competent cells. All transformants from E. coli 
10 ESI 301 mutS were pooled and £. coli SCSI 10 were pooled. A plasmid pool was isolated 
from each pooled library, and this pool was used to transform S. lividans TK23 
protoplasts to form a mutant library for deacylase activity screening. 
Transformants from the S. lividans TK23 libraries were screened for ECB deacylase 
activity with an in situ plate assay. Transformed protoplasts were allowed to regenerate 
1 5 on R2YE agar plates for 24 hr at 30°C and to develop in the presence of thiostrepton for 
48 hours. When the colonies grew to the proper size, 6 ml of 0.7% agarose solution 
containing 0.5 mg/ml ECB in 0.1 M sodium acetate buffer (pH 5.5) was poured on top of 
each R2YE-agar plate and allowed to develop for 18-24 hr at 30°C. Colonies surrounded 
by a clearing zone larger than that of a control colony containing wild-type plasmid 
20 pSHP 150-2*, were selected for further characterization. 

Selected transformants were inoculated into 20 ml medium containing 
thiostrepton and grown aerobically at 30°C for 48 hours, at which point they were 
analyzed for ECB deacylase activity using HPLC. 100 ^il of whole broth was used for a 
reaction at 30 °C for 30 minutes in 0.1 M NaAc buffer (pH 5.5) containing 10% (v/v) 
25 MeOH and 200 }ag/ml of ECB substrate. The reactions were stopped by adding 2.5 

volumes of methanol, and 20 \x\ of each sample were analyzed by HPLC on a 100 x 4.6 
mm polyhydroxyethyl aspartamide column (PolyLC Inc.. Columbia. MD, USA) at room 
temperature using a linear acetonitrile gradient starting with 50:50 of A:B (A = 93% 
acetonitrile. 0.1% phosphoric acid: B = 70% acetonitrile. 0.1% phosphoric acid) and 
30 ending with 30:70 of A:B in 22 mm at a flow rate of 2.2 ml/min. The areas of the ECB 
and ECB nucleus peaks were calculated and subtracted from the areas of the 
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corresponding peaks from a sample culture of S. lividans containing pIJ702* in order to 
estimate the ECB deacylase activity. 

2.0 ml pre-cultures of positive mutants were used to inoculate 50-ml 
medium and allowed to grow at 30°C for 96 hr. The supematants were further 

5 concentrated to 1/30 their original volume using an Amicon filtration unit (Beverly, MA, 
USA) with molecular weight cutoff of 10 kD. The resulting enzyme samples were 
diluted with an equal volume of 50 mM KH 2 P04 (pH 6.0) buffer and were applied to Hi- 
Trap ion exchange column (Pharmacia Biotech, Piscataway, NJ, USA) . The binding 
buffer was 50 mM KH 2 P0 4 (pH 6.0), and the elution buffer was 50 mM KH 2 P0 4 (pH 6.0) 

10 containing 1 .0 M NaCl. A linear gradient from 0 to 1 .0 M NaCl was applied in 8 column 
volumes with a flow rate of 2.7 ml/min. The ECB deacylase fraction eluting at 0.3 M 
NaCl was concentrated and the buffer was exchanged for 50 mM KH 2 P0 4 (pH 6.0) using 
Centricon-10 units. Enzyme purity was verified by SDS-PAGE using Coomassie Blue 
stain, and the concentration was determined using the Bio-Rad Protein Assay Reagent 

15 (Hercules, C A, USA). 

A modified HPLC assay was used to determine the activities of the ECB 
deacylase mutants on ECB substrate (84). Four \ig of each purified ECB deacylase 
mutant was used for activity assay reaction at 30°C for 30 minutes in 0.1 M NaAc buffer 
(pH 5.5) containing 10% (v/v) MeOH and different concentrations of ECB substrate. 
20 Assays were performed in duplicate. The reactions were stopped by adding 2.5 volumes 
of methanol, and the HPLC assays were carried out as described above. The absorbance 
values were recorded, and the initial rates were calculated by least-squares regression of 
the time progress curves from which the Km and the kcat were calculated. 

Activities as a function of pH were measured for the purified ECB 
25 deacylases at 30°C at different pH values: 5, 5.5 and 6 (0.1 M acetate buffer); 7, 7.5, 8 
and 8.5 (0.1 M phosphate buffer); 9 and 10 (0.1 M carbonate buffer) using the HPLC 
assay. Stabilities of purified ECB deacylases were were determined at 30°C in 0.1 M 
NaAc buffer (pH 5.5) containing 10% methanol. Samples were withdrawn at different 
time intervals, and the residual activity was measured in the same buffer with the HPLC 
30 assay described above. 
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B. Results 

Fig. 1 1 shows that after one round of applying this heteroduplex repair 
technique on the mutant M7-2 and Ml 6 genes, one mutant (Ml 5) from about 500 original 
transformants was found to possess 3.1 times the specific activity of wild-type. 
Wild type and evolved M 1 5 ECB deacy lases were purified and their kinetic parameters 
for deacylation of ECB were determined by HPLC. The evolved deacylases Ml 5 has an 
increased catalytic rate constant, by 205%. The catalytic efficiency (kc a i/K m ) of M20 
is enhanced by a factor of 2.9 over the wild-type enzyme. 

Initial rates of deacylation with the wild type and Ml 5 at different pH 
values from 5 to 10 were determined at 200 ^g/ml of ECB. The recombined Ml 5 is 
more active than wild type at pH 5-8. Although the pH dependence of the enzyme 
activity in this assay is not strong, there is a definite shift of 1 .0-1.5 units in the optimum 
to lower pH, as compared to wild type. 

The time courses of deactivation of the purified ECB deacylase mutant 
5 M15 was measured in 0.1 M NaAc (pH 5.5) at 30°C. No significant difference in 
stability was observed between wild type and mutant Ml 5. 

The DNA mutations with respect to the wild type ECB deacylase sequence 
and the positions of the amino acid substitutions in the evolved variants M7-2, Ml 6 and 
Ml 5 are summarized in Figure 12. 

0 The heteroduplex recombination technique can recombine parent 

sequences to create novel progeny. Recombination of the M7-2 and Ml 6 genes yielded 
Ml 5, whose activity is higher than any of its parents (Fid. 13). Of the six base 
substitutions in M15, five (at positions ct50. <xl71. 357, P129 and p340) were inherited 
from M7-2, and the other one ((330) came from Ml 6. 

5 This approach provides an alternative to existing methods of DNA 

recombination and is particularly useful in recombining large genes or entire operons. 
This method can be used to create recombinant proteins to improve their properties or to 
study structure-function relationship. 
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EXAMPLE 3. Novel Thermostable Bacillus Subtilis Subtilisin E Variants 

This example demonstrates the use in vitro heteroduplex formation 
followed by in vivo repair for combining sequence information from two different 
sequences in order to improve the thermostability of Bacillus subtilis subtilisin E. 

Genes RC1 and RC2 encode thermostable B. sublilis subtilisin E variants 
(88). The mutations at base positions 1 107 in RC1 and 995 in RC2 (Figure 14), giving 
rise to amino acid substitutions Asn218/Ser (N218S) and Asnl81/Asp (N181 ID), lead to 
improvements in subtilisin E thermostability; the remaining mutations, both synonymous 
and nonsynonymous, have no detectable effects on thermostability. At 65°C, the single 
variants N181D and N218S have approximately 3 -fold and 2-fold longer half-lives, 
respectively, than wild subtilisin E. and variants containing both mutations have half-lives 
that are 8-fold longer (88). The different half-lives in a population of subtilisin E variants 
can therefore be used to estimate the efficiency by which sequence information is 
combined. In particular, recombination between these two mutations (in the absence of 
5 point mutations affecting thermostability) should generate a library in which 25% of the 
population exhibits the thermos/ability of the double mutant. Similarly, 25% of the 
population should exhibit wild-type like stability, as N181D and N218S are eliminated at 
equal frequency. We used the fractions of the recombined population as a diagnostic 

0 A. Methods 

The strategy underlying this example is shown in Fig. 15. 
Subtilisin E thermostable mutant genes RC1 and RC2 (Fig. 14) are 986-bp 
fragments including 45 nt of subtilisin E prosequence, the entire mature sequence and 1 13 
nt after the stop codon. The genes were cloned between Bam HI and Nde I in £. colilB. 
25 subtilis shuttle vector pBE3, resulting in pBE3-l and pBE3-2, respectively. Plasmid 
DNA pBE3-l and pBE3-2 was isolated from E.coli SCSI 10. 

About 5.0 ^g of ummethylated pBE3-l and pBE3-2 DNA were digested 
with Bam HI and Nde L respectively, at 37°C for 1 hour. After agarose gel separation, 
equimolar concentrations (2.0nM) of the linearized unmethylated pBE3-l and pBE3-2 
30 were mixed in 1 x SSPE buffer ( 1 80 mM NaCL 1 .0 mM EDTA. 1 0 mM NaH 2 P0 4 . pH 

7.4). After heating at 96°C for 10 minutes, the reaction mixture was immediately cooled 
at 0°C for 5 min. The mixture was incubated at 68°C for 2 hr for heteroduplexes to form. 
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One microliter of the reaction mixture was used to transform 50 ul of E. 
coli ES 1301 mutS. E. coli SCSI 10 and E. coli HB101 competent cells. 

The transformation efficiency with E. coli HB101 competent cells was 
about ten times higher than that of E.coli SCSI 10 and 15 times higher than that of E. coli 
5 ES1301 mutS. But in all these cases, the transformation efficiencies were 10-250 times 
lower than that of the transformation with closed, covalent and circular control pUC19 
plasmids. 

Five clones from E. coli SCSI 10 mutant library and five from E. coli 
ESI 301 mutS library were randomly chosen, and plasmid DNA was isolated using a 
1 0 QIAprep spin plasmid miniprep kit for further DNA sequencing analysis. 

About 2.000 random clones from E. coli HB101 mutant library were 
pooled and total plasmid DNA was isolated using a Q1AGEN-100 column. 0.5-4.0 ug of 
the isolated plasmid was used to transform Bacillus subtilis DB428 as described 
previously (88). 

1 5 About 400 transformants from the Bacillus subtilis DB428 library were 

subjected to screening. Screening was performed using the assay described previously 
(88), on succinyl-Ala-Ala-Pro-Phe-p-nitroanilide. B. subtilis DB428 containing the 
plasmid library were grown on LB plates containing kanamycin (20 ug/ml) plates. After 
18 hours at 37°C single colonies were picked into 96-well plates containing 200 ul 
20 SG/kanamycin medium per well. These plates were incubated with shaking at 37°C for 
24 hours to let the cells to grow to saturation. The cells were spun down, and the 
supematants were sampled for the thermostability assay. 

Two replicates of 96-well assay plates were prepared for each growth plate 
by transferring 10 ul of supernatant into the replica plates. The subtilisin activities were 
25 then measured by adding 100 ul of activity assay solution (0.2 mM succinyl-Ala-Ala-Pro- 
Phe-p-nitroanilide. 100 mM Tris-HCl, 10 mM CaCl 2 , pH 8.0. 37°C). Reaction velocities 
were measured at 405 nm to over 1.0 min in aThermoMax microplate reader (Molecular 
Devices. Sunnyvale CAV Activity measured at room temperature was used to calculate 
the fraction of active clones (clones with activity less than 1 0% of that of wild type were 
30 scored as inactive). Initial activity (Ai) was measured after incubating one assay plate at 
65°C for 10 minutes by immediately adding 100 ul of prewarmed (37°C) assay solution 
(0.2 mM succinyl-Ala-Ala-Pro-Phe-p-nitroanilide. 100 mM Tris-HCl. pH 8.0. 10 mM 
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CaCh, pH 8.0) into each well. Residual activity (Ar) was measured after 40 minute 
incubation. 

B. Results 

In vitro heteroduplex formation and in vivo repair was carried out as 
5 described above. Five clones from E. coli SCSI 10 mutant library and five from £. coli 
ES 1301 mutS libraries were selected at random and sequenced. Fig. 14 shows that four 
out of the ten clones were different from the parent genes. The frequency of occurrence 
of a particular point mutation from parent RC1 or RC2 in the resulting genes ranged from 
0% to 50%, and the ten point mutations in the heteroduplex have been repaired without 
10 strong strand-specific preference. 

Since none of the ten mutations locates within the dcm site, the mismatch 
repair appears generally done via the £. coli long-patch mismatch repair systems. The 
system repairs different mismatches in a strand-specific manner using the state of N6- 
methylation of adenine in GATC sequences as the major mechanism for determining the 
15 strand to be repaired. With heteroduplexes methylated at GATC sequences on only one 
DN A strand, repair was shown to be highly biased to the unmethylated strand, with the 
methylated strand serving as the template for correction. If neither strand was 
methylated, mismatch repair occurred, but showed little strand preference (23, 24). These 
results shows that it is preferable to demethylate the DNA to be recombined to promote 
20 efficient and random repair of the heteroduplexes. 

The rates of subtilisin E thermo-inactivation at 65 °C were estimated by 
analyzing the 400 random clones from the Bacillus subtilis DB428 library. The 
thermostabilities obtained from one 96-well plate are shown in Figure 16, plotted in 
descending order. About 12.9% of the clones exhibited thermostability comparable to the 
25 mutant with the N181D and N21 8S double mutations. Since this rate is only half of that 
expected for random recombination of these two markers, it indicates that the two 
mismatches at positions 995 and 1 107 within the heteroduplexes have been repaired with 
lower position randomness. 

Sequence analysis of the clone exhibiting the highest thermostability 
30 among the screened 400 transformants trom the E. coli SCSI 10 heteroduplex library 
confirmed the presence of both N181D and N2 1 8S mutations. Among the 400 
transformants from the B.sublilis DB428 library that were screened, approximately 91% 
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p u„ M on Qia g e„ co.ua™ and Rested w„n toW and M enrvmes. D*«ed 
products were lig a t ed into ,*I and £ coRI Rested pGFP vector. 

dUTP was added into PCR reason at fmal concentrates 200 uM. 40uM. 
8 uM, 1 -6 uM, 0.32 uM. PCR reaction and subsequent ctoatag procedures were 

performed as described above. 

P Recent plasmids were transformed in. XL10 E. coU stram by a 

lifted che mi cal transformation method. Cells were plated on ampicillin containing LB 
modified chemi incuba tion at room temperature or 

agar plates and grown overnight at 37 U toiiow y 

at 4°C until fluorescence developed. 
B. Results . 

i ^^^ a£iigglioJ}mie ^ h i n mion efficienc y,. 
Two experiments have been performed to test the effect of breaks in the 
, , n the efficiency of recombination. In one experiment heteroduplex 
DNA heteroduplex on the efficiency o single . stran d breaks and was 

5 piasmid was treated with DNA ligase to close all existing s g 

transformed in identical conditions as an unligated sample (see Table 1). hp 
TpL show up to 7-fold improvement in recombination efficiency over unligated 

SamPleS ' in another experiment, dUTP was added into PCR reaction to introduce 

«r hv uracvl N-elvcosvlase in the host cells. Table 
20 additional breaks into DNA upon repair b ^ ^ V ; ^ ^ of 

■> shows that dUMP incorporation sigmficantly suppressed 
suppression increasing with increased dUTP concentration. 

~ rfTrrt nf p 1 — ;H ^ nn thP pfflciencY ° f 
-,c hrttrnduple * formation. 

P,asmid size was a significant factor affecting recombination efficency. 
, 'A „GFP (3 3 kt» and a Bacillus shuttie vector pCT. (about 9 Kb) were used 

40% of circular plasmid. the shuttle vector yi 
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WHAT IS CLAIMED IS: 



1 1 . A method for evolving a polynucleotide toward acquisition of a 

2 desired property, comprising 

3 (a) incubating a population of parental polynucleotide variants under 

4 conditions to generate annealed polynucleotides comprising heteroduplexes; 

5 (b) exposing the heteroduplexes to a cellular DNA repair system to 

6 convert the heteroduplexes to parental polynucleotide variants or recombined 

7 polynucleotide variants; 

8 (c) screening or selecting the recombined polynucleotide variants for the 

9 desired property. 

1 2. The method of claim 1 , wherein the heteroduplexes are exposed to 

2 the cellular DNA repair system in vitro. 

1 3 . The method of claim 2, wherein the cellular DNA repair system 

2 comprises cellular extracts. 

1 4. The method of claim 1 , further comprising introducing the 

2 heteroduplexes into cells, whereby the heteroduplexes are exposed to the DNA repair 

3 system of the cells in vivo. 

1 5. The method of claim 4, wherein the annealed polynucleotides 

2 further comprise homoduplexes and the introducing step selects for transformed cells 

3 comprising the heteroduplexes relative to transformed cells comprising homoduplexes. 

1 6. The method of claim 4, wherein a first polynucleotide variant is 



2 provided as a component of a first vector, and a second polynucleotide variant is provided 

3 as a component of a second vector, and the method further comprises converting the first 

4 and second vectors to linearized forms in which the first and second polynucleotide 

5 variants occur at opposite ends, whereby in the incubating step single-stranded forms of 

6 the first linearized vector reanneal with each other to form linear first vector, single- 

7 stranded forms of the second linearized vector reanneal with each other to form linear 

8 second vector, and single-stranded linearized forms of the first and second vectors anneal 

9 with each to form a circular heteroduplex bearing a nick in each strand, and the 
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intr odu Cl ng step selects for transformed cells composing the circular 
heteroduplexes relative to the linear first and second vector. 

7 . Tne method of claim 6. wherein the first and second vectors are 
converted to linearized forms by PCR. 

8 The metftod of data, 6, wherem .he firs, and second vecors are 
convened to Unearned forms by d.gestmn wiih first and second resuiction enzyme, 

9 The method of claim U wherein the populat.cn of polynucleotide 

a a f.rm and the method further comprising 

the annealing step. 

10 The method of claim 1, wherein the convening step comprises: 
inducting asymmetric amplification of the firs, and second double 
• ♦ omnlifv a first strand of the first polynucleotide 

second strands anneal in me incubating step to form a heteroduple*. 

, ! The method of claim 10. where.n the firs, and second double- 
strand ed polvncueltoides vanants are provided ,n vector-free form, and the method 
funher comprises incorporating the heteroduplex into a vector. 

12 The method of claim 4 wherein the population of polynucleotides 

components of firs, and s ^ ^ ^ ^ 
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strand, and the introducing step selects for transformed cells comprises the 

2 circular heteroduplexes relative to the linear first and second vector. 

1 13. The method of claim 4, further comprising sealing nicks in the 

2 heteroduplexes to form covalently-closed circular heteroduplexes before the introducing 

3 step. 

1 14. The method of claim 1 1 , wherein the first and second 

2 polynucleotides are obained from chromosomal DNA.. 

1 15. The method of claim 1 , further comprising repeating steps (a)-(c) 

2 whereby the incubating step in a subsequent cycle is performed on recombinant variants 

3 from a previous cycle. 

1 16. The method of claim 1 , wherein the polynucleotide variants encode 

2 a polypeptide. 

1 17. The method of claim 1 , wherein the population of polynucleotide 

2 variants comprises at least 20 variants. 

1 18. The method of claim 1 , wherein the population of polynucleotide 

2 variants are at least 10 kb in length. 

1 19. The method of claim L wherein the population of polynucleotide 

2 variants comprises natural variants. 

1 20. The method of claim 1 , wherein the population of polynucleotides 

2 comprises variants generated by mutagenic PCR. 

1 21 . The method of claim 1 , wherein the population of polynucleotide 

2 variants comrises variants generated by site directed mutagenesis. 

1 22. The method of claim 1 . wherein the cells are bacterial cells. 

1 23 . The method of claim 1 . further comprising at least partially 

2 demethylating the population of variant polynucleotides. 
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1 24. The method of claim 23. whether the at least partially 

2 demethylating step is performed by PCR amplification of the population of variant 

3 polynucleotides. 

1 25. The method of claim 23. wherein the at least partially 

2 demethylating step is performed by amplification of the population of variant 

3 polynucelotides in host cells. 

1 26. The method of claim 25, wherein the host cells are defective in a 

2 gene encoding a methylase enzyme. 

1 27. The method of claim 1, wherein the population of variant 

2 polynucleotide variants comprises at least 5 polynucleotides having at least 90% sequence 

3 identity with one another. 

1 28. The method of claim 1 , further comprising isolating a screened 

2 recombinant variant. 

1 29. The method of claim 28, further comprising expressing a screened 

2 recombinant variant to produce a recombinant protein. 

1 30. The method of claim 29 further comprising formulating the 

2 recombinant protein with a carrier to form a pharmaceutical composition. 

1 31. The method of claim 1 , wherein the polynucleotide variants encode 

2 enzymes selected from the group consisting of proteases, lipases, amylases, cutinases, 

3 cellulases, amylases, oxidases, peroxidases and phytases. 

1 32. The method of claim 1, wherein the polynucleotide variants encode 

2 a polypeptide selected from the group consisting of insulin. ACTH. glucagon. 

3 somatostatin, somatotropin, thymosin, parathyroid hormone, pigmentary hormones. 

4 somatomedin, erthropoietin. luteinizing hormone, chorionic gonadotropin, hyperthalmic 

5 releasing factors, antidiuretic hormones, thyroid stimulating hormone, relaxin. interferon. 

6 thrombopoietin (TPCM. and prolactin. 
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1 33. The method of claim 1, wherein the polynucleotide variants 

2 encode a plurality of enzymes forming a metabolic pathway. 

1 34. The method of claim L wherein the polynucleotide variants are in 

2 concatemeric form. 
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, CTCC^^CCCCAGCTGTTCCTOG^^TCCCGGCCGCCC^GCCGCCOTCCCaOTC 

^ , OCCGCCGCGGCGCCGATCG^TCGXCCCTrCGrCGTGCCGCAG^CGCCCTGCCGCTC 

121 TGCGGCGGCAGCCGGCCGCCCCTGCTCGCCTCGGCGATGCTCGGCGCGCTGCTGGTGGTC 

1C1 ggCGCCGACCTGGTCGCTCAGATCGTGGTGGCGCCGAAGGAGCTGCCGGTCGGCCTGCTC 

241 ACCGCGATGATCGGCACCCCGTACCTGCTCTGGCTCCTGCTTCGGCGATCAAGAAAGGTG 

„1 ACGGGTGCGGCAAGTCGACGCTGCTCAAGGCGCTCGGCCGGCTGCTGCGCCCGACCGGCG 

4.1 GGCAGGTGCTGCTGGACGGCCGCCGCATCGACCGGACCCCCACCCGTGACGTGGCCCGGG 

541 TCCTCGGCGTGCTGCCGCAGTCGCCCACCGCGCCCGAAGGGCTCACCGTCGCCGACC^G 

6 01 XGATGCGCGGCCGGCACCCGCACCAGACCTGGTTCCGGCAGTGGTCGCGCGACGACGAGG 

721 ACGCCCTCTCCGGCGGTCAGCGCCAGCGCGCCTGGATCAGCATGGCGCTGGCCCAGGGCA 
781 CCGACCTGCTGCTGCTGGACGAGCCGACCACCTTCCTCGACCTGGCCCACCAGATCGACG 
.41 TGCTGGACCTGGTCCGCCGGCTGCACGCCGAGATGGGCCGGACCGTGGTGATGGTGCTGC 

96 1 TCGTGGCGAGCGGGGCGCCGGACGAGGTGCTCACCCCGGCGCTGCTGGAGTCGGTCTTCG 

10B1 CGCGCCCCGCCACCTCGGTGCGGGCCTGAAATCGATGAGCGTGGTTGCTTCATCGGCCTG 
H41 CCGAGCGATGAGAGTATGTGGGCGGTAGAGCGAGTCTCGAGGGGGAGATGCCGCCGTGAC 



:tgaaagcagcagcgatcgccttcggtgtgatcgtggcgaccgc 



1201 GTCCTCGTACATGCGCCT 

K A A A I 

VTGACGGCGGCTATGCGGCCCTGATCCG 



S SYMRLKAAAIA^ GVJ " 



GTGCCGTCACCCGCTTCCGGCAGGGAACA1 



1261 AGCO 

~ crREHDGGYAALXR 
23 avPSPASGRE"" 



1321 CCGGGCCTCGTACGGCGTCCCGCACATCACCGCCGACGACTTCGGGAGCCTCGG' 
43 RASYGVP 



TTTCGG 

KITADDF6SLGFG 



\ 3 

1441 CAACGGTGAGCGGTCGCGGTGGTTCGGXGCG.CCGGGCCGGACGACGCCGATGTGCGCAG 

Fie. 13A 



WO 99/29902 



15 / 21 



PCT/US98/25698 



83 NGERSRWFGATGPDDADVRS 

15 01 CGACCTCTTCCACCGCAAGGCGATCGACGACCGCGTCGCCGAGCGGCTCCTCGAAGGGCC 
103 DLFHRKA I DDRVAERLLEGP 

1561 CCGCGACGGCGTGCGGGCGCCGTCGGACGACGTCCGGGACCAGATGCGCGGCTTCGTCGC 
123 K D G V R A P SDDVRDQMRGFVA 

1621 CGGCTACAACCACTTCCTACGCCGCACCGGCGTGCACCGCCTGACCGACCCGGCGTGCCG 
143 GYNHFLRRTGVHRLTDPACR 

1681 CGGCAAGGCCTGGGTGCGCCCGCTCTCCGAGATCGATCTCTGGCGTACGTCGTGGGACAG 
163 GKAWVRPLSEIDLWRTSWDS 

17 41 CATGGTCCGGGCCGGTTCCGGGGCGCTGCTCGACGGCATCGTCGCCGCGACGCCACCTAC 
183 MVRAG SGALLDGIVAATP PT 

18 01 AGCCGCCGGGCCCGCGTCAGCCCCGGAGGCACCCGACGCCGCCGCGATCGCCGCCGCCCT 
203 AAGPASAPEAPDAAAIAAAL 

18 61 CGACGGGACGAGCGCGGGC ATCGGCAGCAACGCGT ACGGCCTCGGCGCGCAGGCCACCGT 
223 DGTSAGI GSNAYGLGAQATV 

1921 GAACGGCAGCGGGATGGTGCTGGCCAACCCGCACTTCCCGTGGCAGGGCGCCGCACGCTT 
243 NGSGKVLANPHFPWQGAARF 

1981 CTACCGGATGC ACCTCAAGGTGCCCGGCCGCTACGACGTCGAGGGCGCGGCGCTGATCGG 
263 YRMHLKVPGRYDVEGAALIG 

2 041 CGACCCGATCATCGGGATCGGGCACAACCGCACGGTCGCCTGGAGCCACACCGTCTCCAC 
2B3 DP I 2 G I GHNRTVAWSHTVST 

2101 CGCCCGCCGGTTCGTGTGGC ACCGCCTGAGCCTCGTGCCCGGCGACCCCACCTCCTATTA 

3 03 ARRFVWHRLSLVPGDPTSYY 
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1 6 1 CCTCGACGGCCGGCCCGAGCGGATGCGCGCCCGCACGGTCACGGTCCAGACCGGCAOCGG 



2 

323 



VDGRPERMRARTVTVQTGSG 



2221 CCCGGTCAGCCGCACCTTCCACGACACCCGCTACGGCCCGGTGGCCGTGATGCCGGGCAC 
343 PVSRTTHDTRYGPVAVMPGT 

2261 CTTCGACTGGACGCCGGCCACCGCGTACGCCATCACCGACGTC AACGCGGGCAAC AACCG 
363 FDWTPATAVAITDVNAGNNR 

2 3 4 1 CGCCTTCGACGGOTGGCTGCGGATGGGCCAGGCCAAGGACGTCCGGGCGCTCAAGGCGGT 
383 XFDO WLRKSQAK''"dVRALKAV 

2401 CCTCGACCGGCACCAGTTCCTGCCCTGGGTCAACGTGATCGCCGCCGACGCGCGGGGCGA 
403 LDRHQFLPWVNVIAADARGE 

2461 GGCCCTCTACGGCGATCATTCGGTCGTCCCCCGGGTGACCGGCGCGCTCGCTGCCGCCTG 
423 ALYGDHSVVPRVTGALAAAC 

CATCCCGGCGCCGTTCCAGCCGCTCTACGCCTCCAGCGGCCAGGCGGTCCTGGACGGTTC 



2521 
443 I P A 



PFQPLYASSGQAVLDGS 



CGCCGACCCCGACGCCGCGGTCCCGGGCATTCTCGGCCC 
DCALGADPDAAVPGILGP 



2581 CCGGTCGGACTGCGCGCTCG3 1 
4 63 R S 



2641 GGCGAGCCTGCCGGTGCGGTTCCGCGACGACTACGTCACCAACTCCAACGACAGTCACTG 
483 A S L P V R F R D D Y V T H S N D S H H 

2701 GCTGGCCAGCCCGGCCGCCCCGCTGGAAGGCTTCCCGCGGATCCTCGGCAACGAACGCAC 
503 I.ASPAAPLEGFPRILGNERT 

27 61 CCCGCGCAGCCTGCGCACCCGGCTCGGGCTGGACCAGATCCAGCAGCGCCTCGCCGGCAC 
S23 pRSLRTRLGLDQlQQR^AGT 



1 GGACGGTCTGCCCGGCAAGGGCTTCACCACCGCCCGGCTCTGGCAGGTCATGTTCGOCAA 
S43 D G L P G K G F T T A R X, W Q V K F G N 



282 



Fig. 13C 
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28el CCOGATG^CGGCGCCGAACTCGCCCGCGACGACCT^TCGCGCTCTGCCGCCCCCAGCC 
5,3 R H H G A E L A R D D L V A t C R R Q * 

2941 GACCGCGACCGCCTCGAACGGCGCGATCGTCGACCTCACCGCGGCCTGCACGGCGCTGTC 
SS 3 T A Z A S H G A I V D L T A A C T A L S 
3001 CCGCTTCGATGAGCGTGCCGACCTGGACAGCCGGGGCGCGCACCTGTTCACCGAGTTCGC 



603 R F D 



E R A D L D S R Gr .A H L F T E F A 



1 CCTCGCGGGCGGAATCAGGTTCGCCGACACCTTCGAGGTGACCGATCCGGTACGCACCCC 



306 

623 L A G G I 



RFADTFE VTDPVRTP 



PRVRTALADAVQR 



643 R R L N T T D 



31.X GCTCGCCGGCATCCCCCTCGACGCGAAGCTGGGAGACATCCACACCGACAGCCGCGGCGA 
663 L A G I PL 



DAKLGDIHTDSRGE 



3241 ACGGCGCA' 
6B3 R R I P I 



TCCCCATCCACGGTGGCCGCGGGGAAGCAGGCACCTTCAACGTGATCACCAA 



HGGRGEAGTFNVITK 



3301 CCCGCTCGTGCCGGGCGTGGGATACCCGCAGGTCGTCCACGOAACATCGTTCGTGATGGC 
703 P L V P G 



VGYPQVVHGTSFVMA 



3361 CGTCGAACTCGGCCCGCACGGCCCGTCGGGACGGCAGATCCTCACCTATGCGCAGTCGAC 
723 v E I. 0 P « G P S G R Q X L T V A Q S T 
3421 GAACCCGAACTCACCCTGGTACGCCGACCAGACCGTGCTCTACTCGCGGAAGGGCTGGGA 
743 N P N S P W V A D Q T V L Y S R K G « D 

763 T X K Y T E A Q I A - O P K , R V Y R V 



3 541 GGCACAGCGGGCA! 



:G CTGACCCACGTCACGCCGGCTCGGCCCGTGCGGGGGCGCAGGGCG 
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793 A Q R G R 

3601 ccGATCGTCTCTGCATCGCCGGTCAGCCGGGGCCTGCGTCGACCGGCGGCGGCCGGTCGA 

3 G 6 1 CGCCCGCGTCCCGGCGCAGCGACTGGCTGAAGCGCCAGGCGTCGGCGGCCCGGGGCAGGT 

3721 TGTTGAACATCACGTACGCCGGGCCGCCGTCGAGGATGCCGGCGAGGTGTGCCAGCTCGG 

3 7 81 C ATCCGTGTACACATGCCGGGCGCCGGTGATGCCGTGCAGCCGGTAATAGGCCATCGGCG 

3 841 TCAGACTGCGGCGCAGGAACGGGTCGGCGGCGTGGGTCAGGTCCAGCTCCTGGCACAAGC 

3 901 CCTCGACCACCTCGTCCGGCCACGGGCCGCGCGGCTCCCACAACAGCCGGACACCGGCCG 
3951 GCCGGCGCGCTCGGGCGCAGAACTCACGCAGTCGCGCGATGGCGGGTTCGGTCGGCCGGA 

4 021 AACTCGCCGGGCACTGCAG 



Fig. 13E 
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